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Abstract 



We report on the observation of i? ^ D*it^-k vr vr'' decays. The branching 
ratios for D*+ and D*^ are (1.72±0.14±0.24)% and (1.80±0.24±0.27)%, re- 
spectively. Each final state has a 0*0071" component, with branching ratios 
(0.29±0.03±0.04)% and (0.45±0.10±0.07)% for the D*+ and D*^ modes, re- 
spectively. We also observe B Du;7r~ decays. The branching ratios for 
and Z)° are (0.28±0.05±0.04)% and (0.41±0.07±0.06)%, respectively. A spin 
parity analysis of the ujtt" system in the Dujtt" final state shows a preference 
for a wide 1~ resonance. A fit to the lvtt~ mass spectrum finds a central mass 
of (1349±25t5°) MeV and width of (547±86t4^) MeV. We identify this object 
as the /9(1450) or the p' . 
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I. INTRODUCTION 



Currently most B meson decays are unknown. The decay width is comprised of hadronic 
decays, leptonic decays and semileptonic decays. Leptonic decays are predicted to be very 
small, ~ 10~^ in branching fraction. The semileptonic branching ratio for B Xe~v, 
X/i~z/, and Xt~v totals approximately 25% [Q. The remainder must come from hadronic 
decays, where the measured exclusive branching ratios total only a small fraction of the 
hadronic width. 

To be more explicit, the measured hadronic decay modes for the B including D^irm^), 
D*+{n7r-), where 3 > n > 1, D+^*W~'-*\ and J/ip exclusive totals only about 10% The 
B^ modes total about 12%. 

Yet, understanding hadronic decays of the B is crucial to ensuring that decay modes 
used for measurement of CP violation truly exhibit the underlying quark decay mechanisms 
expected theoretically. 

It is also interesting to note that the average charged multiplicity in a hadronic B^ decay 
is 5.8±0.1 0. Since this multiplicity contains contributions from the decay of or D*~^ 
normally present in 5^ decay, we expect a sizeable, approximately several percent, decay rate 
into final states with four pions The seen D^*\mi)~ final states for n < 3 are consistent 
with being quasi-two-body final states. For n of two the p~ dominates, while for n of three 
the Qi dominates 0. These decays appear to occur from a simple spectator mechanism 
where the virtual W~ materializes as a single hadron: vr^, or aj~. The decay rates can be 
understood in a simple "factorization" model where the decay rate is given by the product 
between two currents, one between the B and the D and the other given by the virtual W~ 
transforming into the light hadron of interest 

In this paper we investigate final states for n of 4. We will show a large signal for 
the D*^ 71^ 71^ 71^ 71^ final state in Section |T|. In Section |^ we will show that a substantial 
fraction, ~20%, arise from D*^uj7c^ decays and that the uJ7i~ mass distribution has a resonant 
structure around 1.4 GeV with a width of about 0.5 GeV. In Section 0, similar conclusions 
are drawn about the D*^7i~^7[~'7[~7i'^ final state. The same structure is shown to exist in 
Du!7r~ final states (Section |VI1]) and we will use these events to show in Section |VIIF| that 
the spin-parity is most likely 1~. This state is identified as the p', sometimes called the 
p(1450). We perform a generalized Breit-Wigner fit to get the best values of the mass and 
width, discussed in Section |X| and the Appendix. 

Other resonant substructure is searched for, but not found (Section Finally we 
summarize our findings and compare with the predictions of factorization and other models 
in Section [XI| . 

The data sample consists of 9.0 fb~^ of integrated luminosity taken with the CLEO II 
and II. V detectors 0] using the CESR e'^e~ storage ring on the peak of the T(45') resonance 
and 4.4 fb~^ in the continuum at 60 MeV less center-of-mass energy. The sample contains 
19.4 million B mesons. 
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II. COMMON SELECTION CRITERIA 



Hadronic events are selected by requiring a minimum of five charged tracks, total visible 
energy greater than 15% of the center-of-mass energy, and a charged track vertex consistent 
with the nominal interaction point. To reject continuum we require that the Fox- Wolfram 
moment R2 be less than 0.3 [0. 

Track candidates are required to pass through a common spatial point defined by the 
origin of all tracks. Tracks with momentum below 900 MeV/c are required to have an ion- 
ization loss in the drift chamber within 3cr of that expected for their assigned mass]^ (These 
requirements are not imposed on slow charged pions from D*^ decay.) Photon candidates 
are required to be in the "good barrel region," within 45° of the plane perpendicular to 
the beam line that passes through the interaction point, and have an energy distribution in 
the Csl calorimeter consistent with that of an electromagnetic shower. To select vr^'s, we 
require that the diphoton invariant mass be between -3.0 to +2.ba of the 7r° mass, where a 
varies with momentum and has an average value of approximately 5.5 MeV. The two-photon 
candidates are then kinematically fit by constraining their invariant mass be equal to the 
nominal 7r° mass. 

We select and candidates via the decay modes shown in Table |. We require that 
the invariant mass of the D candidates lie within ±2.5(j of the known D masses. The cr's 
are also listed in Table |. The widths vary with the momentum, p, (units of MeV), 
while the widths are not momentum dependent. 

We select D*^ candidates by imposing the addition requirement that the mass difference 
between tt^D^ and combinations is within ±2.5(T of the known mass difference. For the 
we use the same requirement for the ti^D^ decay. The mass difference resolutions are 
0.63 and 0.90 MeV, for the ti^D^ and -n^D^ modes, respectively p. 



TABLE I. Mass Resolutions (cr) in MeV [p in units of MeV) 





K vr+vr^ 


K-TT+ 


DO ^ K-tt+ttO 


D° ^ ET-vr+Tr+vr- 




6.0 


px0.93xl0~3+6.0 


px0.68xl0"2+11.6 


pxO. 92x10^^+4.7 



III. OBSERVATION OF D*+7r+7r vr TfO DECAYS 

A. B Candidate Selection 

We start by investigating the D*+(47r)~ final state.0 The D*^ candidates are pooled with 
all combinations of n^Tf^n^'n-^ mesons. 



Here and throughout this paper a indicates an r.m.s. error. 
^In this paper (47r)^ will always denote the specific combination ir^ Tr~ tt^ . 
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Next, we calculate the difference between the beam energy, Eheam, and the measured 
energy of the five particles, AE. The "beam constrained" invariant mass of the B candidates, 
Mb, is computed from the formula 
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To further reduce backgrounds we define 



xl 




(2) 



where AM^* is the computed D* — mass difference minus the nominal value, AMd is 
the invariant candidate mass minus the known mass and AM^o is the measured 77 
invariant mass minus the known ir^ mass. All vr^'s in the final state are included in the sum. 
The (t's are the measurement errors. We select candidate events in each mode requiring that 
Xb < Cni where Cn varies for each decay decay mode. For the Kmr decay modes we use 
Cn= 12, 8, and 6, respectively. 
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FIG. 1. The B candidate mass spectra for the final state Z)*^7r^7r vr vr" 



-° with K-TT+ (a) 

for AE sidebands and (b) for AE consistent with zero. The curve in (a) is a fit to the background 
distribution described in the text, while in (b) the shape from (a) is used with the normalization 
allowed to float and a signal Gaussian of width 2.7 MeV is added. 



B. Branching Fraction and (47r) Mass Spectrum 

We start with the K^tt^ decay mode. We show the candidate B mass distribution, 

Mb, for AE in the side-bands from -5.0 to -3. Oct and 3.0 to 5.0a on Fig. |l|(a). The AE 
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resolution is 18 MeV (cr). The sidebands give a good representation of the background in the 
signal region. We fit this distribution with a shape given as 



back{r) = pirVl — T^e '^^^^ , 
where r = Mb/5.2895 GeV, and the Pi are parameters given by the fit. 



(3) 



600 



400 



>200 



CM 

c 

iiieoo 



400 



200 



_ 1 1 1 1 1 1 1 1 1 1 
-(a) 


1 1 1 1 1 1 ^ 1 ■ 


— \ — \ — 1 — 1 — 1 — \ — \ — 1 — 1 — 1 — h- 

(b) 








" , , , 1 , , , 1 , , 





600 



400 



>200 



CM 

c 

ii;6oo 



400 



200 



(c) 



—I — I — \ — I — I — I- 
(d) 




H 1 1 — 1 1- 




5.20 5.22 5.24 5.26 5.28 5.30 5.20 5.22 5.24 5.26 5.28 5.30 

Mg (GeV) Mg (GeV) 

FIG. 2. The B candidate mass spectra for the final state D*^ tt^ tt^ ir^ tt^ , the left-side plots are 
for K~Tr~^7r^ (a) AE sidebands, (b) for AE consistent with zero; the right-side plots are for 

'Tr'^TT^ir (c) AE sidebands, (d) for AE consistent with zero. The curves in the top plots, 

(a) and (c), are fits to the background distribution described in the text, while in the bottom plots, 

(b) and (d), the shapes from (a) and (c) are used with the normalization allowed to float and a 
signal Gaussian of width 2.7 MeV is added. 

We next view the Mb distribution for events having AE within 2a around zero in 
Fig. |I](b). This distribution is fit with a Gaussian Signal function of width 2.7 MeV and 
the background function found above whose normalization is allowed to vary. The Gaussian 
signal width is found from Monte Carlo simulation. The largest and dominant component 
results from the energy spread of the beam. We find 358±29 events in the signal peak. 

We repeat this procedure for the other two Z)° decay modes. The Mb spectrum for the 
AE sidebands and the signal region is shown in Fig. 0. The AE resolution is 22 MeV in 



the K- 



mode and 18 MeV in the K tt+tt+tt mode. Signal to background ratios are 



worse in these two modes, but the significance in both modes is quite large. The numbers 
of signal events are shown in Table |I[ 

We choose to determine the branching fraction using only the K^tt^ decay mode 

because of the relatively large backgrounds in the other modes and the decreased systematic 
error due to having fewer particles in the final state. In order to find the branching ratio 
we use the Monte Carlo generated efficiency, shown in Fig. ^ as a function of (47r)~ mass. 
The efficiency falls off at larger (47r)~ masses because the detection of the slow vr"*" from the 
D*~^ decay becomes increasingly more difficult. Since the efficiency varies with mass we need 
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TABLE II. Event numbers for the D*+7r+7r vr Tr*^ final state 



Decay Mode 


Fitted # of events 




358±29 




543±49 




329±41 



to determine the (47r)~ mass spectrum. To rid ourselves of the problem of the background 
shape, we fit the B candidate mass spectrum in 50 MeV bins of (47r)~ mass. (The mass 
resolution is approximately 12 MeV.) The resulting (47r)~ mass spectrum is shown in Fig. ^ 
There are indications of a low- mass structure around 1.4 GeV, that will be investigated 
further in this paper. 

1630301-002 
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FIG. 3. The efficiency for the final state D*+ n+ -k' -k^ , with K--K+. 
We find 

^ D*+7r+7r+7r-7r°) = (1.72 ± 0.14 ± 0.24)% . (4) 

The systematic error arises mainly from our lack of knowledge about the tracking and vr'^ 
efficiencies. We assign errors of ±2.2% on the efficiency of each charged track, ±5% for the 
slow pion from the D*"*", and ±5.4% for the 7r°. The error due to the background shape is 
evaluated in three ways. First of all, we change the background shape by varying the fitted 
parameters by la. This results in a change of ±3%. Secondly, we allow the shape, p2, to 
vary (the normalization, pi, was already allowed to vary). This results in 3.8% increase in 
the number of events. Finally, we choose a different background function 

hack'{r) = pirVl — (1 + + p^r"^ + P4r^] , (5) 
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FIG. 4. The invariant mass spectra of vr+vr^Tr^vr'^ for the final state D*^ tt^ tt^ tt' tt^ , with 
Z?'' ^ K^TT^, found by fitting the B yield in bins of Att mass. 



and repeat the fitting procedure. This results in a 3.7% decrease in the number of events. 
Taking a conservative estimate of the systematic error due to the background shape we 
arrive at ±3.8%. We use the current particle data group values for the relevant D*^ and 
D° branching ratios of (68.3±1.4)% {D*+ 7r+L)°) and (3.85±0.09)% (D° K-7v+), 
respectively The relative errors, 2.0% for the D*^ branching ratio and 2.3% for the are 
added in quadrature to the background shape error, the 7r° detection efficiency uncertainty 
and the tracking error. The total tracking error is found by adding the error in the charged 
particle track finding efficiency linearly for the 5 "fast" charged tracks and then in quadrature 
with the slow pion from the D*~^ decay. The total systematic error is 14%. 

We wish to search for narrow structures. However, we cannot fit the B mass spectrum 
in small (47r)~ mass intervals due to a lack of statistics. Thus we plot the (47r)~ mass 
for events in the Mb peak for the K'n^ mode and the sum of all three modes in 

Fig. 1^. We also plot two background samples: events at lower Mb (5.203 - 5.257 GeV) 
and those in the /S.E sideband separately. First we view the plots in the canonical 50 MeV 
bins. Both background distributions give a consistent if somewhat different estimates of 
the background shape. (Each background distribution has been normalized to the absolute 
number of background events as determined by the fit to the Mb distribution.) In any case 
no prominent narrow structures appear in the histograms for the 10 MeV binning. 

The most accurate distribution of 47r~ mass is obtained by using the data in all three 
decay modes. The 47r~ mass distribution shown in Fig. |^ was found by fitting the Mb 
candidate mass distributions summed together. The distribution has been corrected for 
efficiency as a function of mass. There is an additional 14% systematic scale uncertainty on 
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FIG. 5. The invariant mass spectra of Tr'^Tr^vr^vr'^ for the final state i5*^7r^7r^7r^7r'^, with 
—>■ K^TT^ (upper left), and the sum of all three decay modes (upper right). Events are 
selected by being within 2a of the B mass. The solid histogram is the background estimate from 
the Mb lower sideband and the dashed histogram is from the AE sidebands; both are normalized 
to the fitted number of background events. The same distributions in smaller bins (lower plots). 
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FIG. 6. The efficiency-corrected background-removed invariant mass spectra of vr^vr^Tr^vr'^ for 
the final state D*^ tt^ tt^ ir^ , for the sum of ah three decay modes. (There is an additional 
scale uncertainty of 14%.) 



IV. THE B D*+ujTT- DECAY 

To investigate the composition of the (47r)~ final state, we now investigate the tt+tt^tt'^ 
mass spectrum for the events in the B peak. All three decay modes are used. We show 
the 7r"'"7r~7r° invariant mass distribution for events in the B mass peak in Fig. ^ (there are two 
combinations per event). A clear signal is visible at the u. The histograms on the figure are 
for events either in the lower Mb range, from 5.203 GeV to 5.257 GeV, or in the previously 
defined AE sidebands; no u signal is visible. 

The purity of the u sample can be further improved by restricting candidates to certain 
regions of the Dalitz plot of the decay products. We define a cut on the Dalitz plot as follows. 
Let To, T+ and T_ be the kinetic energies of the pions, and Q be the difference between the 
Lj mass, (equal to 782 MeV), and the mass of the 3 pions. We define two orthogonal 
coordinates X and Y, where 

X = 3To/Q - 1 (6) 
Y = V3{T+-T^)/Q . (7) 

The kinematic limit that defines the Dalitz plot boundary is defined as 

'^boundary ~ '^{-^boundary + ^){Xboundary + 1 + Q.)(l + b/ (^X boundary + 1 — c)) (8) 

where a = Guiq/Q, b = Qrn? / (M^Q), c = 3{M^ — mo)^/ (2M^(5), m is the mass of a charged 
pion and rriQ the mass of the neutral pion. 
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FIG. 7. The invariant mass spectra of vr+vr^vr^ for the final state D*^7r^7r^7r^7r'^ for all three 
decay modes. The solid histogram is the background estimate from the Mb lower sideband 

and the dashed histogram is from the AE sidebands; both are normalized to the fitted number of 

background events. (There are two mass combinations per event.) 



For any set of three pion kinetic energies, we define a variable r, properly scaled to the 
kinematic limit as 



\ -I- 

\ boundary ' boundary 



where the boundary values are found by following the radial vector from (0,0) through (X, Y). 

For events in the B mass peak we show in Fig. | the tt+tt^tt^ invariant mass for three 
different cuts on r. The u signal is purified by restricting r, since the Dalitz plot density for 
a 1~ system decaying into tt^tt^tt^ peaks at r equals zero 0. 

For further analysis we select uj candidates within the 7r+7r^7r° mass window of 782±20 
MeV with r < 0.7. We abandon the xl cut here, as background is less of a problem. In 
Fig. ^ we show the B candidate mass distribution for the D*^uj7t^ final state summing over 
all three D° decay modes. (The signal is fit with the same prescription as before.) There 
are 136±15 events in the peak. 



In Fig. |T0| we show the u!ti~ mass spectrum in the left-side plot. The solid histogram 
shows events from the lower Mb sideband region suitably normalized. The dotted histogram 
shows the background estimate from the AE sidebands, again normalized. In the signal 
distribution there is a wide structure around 1.4 GeV, that is inconsistent with background. 
We re-determine the tuvr" mass distribution by fitting the Mb distribution in bins of uJ7i~ 
mass, and this is shown on the right-side. 

Knowing the cutt" mass dependence of the efficiency we evaluate the branching fraction: 
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FIG. 8. The invariant mass spectra of vr'^vr^vr for the final state D*^ ir^ tt" tt" tt for ah three 
decay modes for three selections on r less than: (a) 1, (b) 0.7 and (c) 0.5. The solid histogram 
is the background estimate from the Mb lower sideband and the dashed histogram is from the AE 
sidebands; both are normalized to the fitted number of background events. 

^ D*+cj7r") = (0.29 ±0.03 ±0.04)% . (10) 

We provisionally label the state at 1432 MeV the A~ and investigate its properties later. 
The U7r~ comprises about 17% of the (47r)~ final state. All of the uJ7i~ final state is consistent 
with coming from A~ decay. 

V. OBSERVATION OF D*^t:+i:--k-tt^ 

We proceed in the same manner as for the reaction with the exception that we use 
the D*^ TT^D^ decay mode and restrict ourselves to the —>■ K^tt^ decay mode only 
due to large backgrounds in the other modes. The is calculated according to equation ^ 
and we use a cut value of 8. The Mb distributions for AE sidebands and signal data are 
shown in Fig. |TT| for the —>■ K~tt^ decay mode. The AE resolution is 18 MeV. We see a 
signal of 195±26 events yielding a branching fraction of 

B{B- D*\+TT-TT-n^) = (1.80 ± 0.24 ± 0.27)% . (11) 
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FIG. 9. The Mb spectra for L'*+u;7r~ for all three decay modes, (a) AE sidebands and 
(b) AE around zero. 



The vr^TT^TT^ mass spectrum shown in Fig. |T2| shows the presence of an u. Selecting on 
the presence of an u with r < 0.7 we show the sideband and signal plots in Fig. 0. (Here 
we do not use the previously defined cut.) The branching ratio, based on 26±6 events is 



B{B- D*"lutt-) = (0.45 ± 0.10 ± 0.07)% 



(12) 



In Fig. |14|we show the uj7i~ mass spectrum. We see an enhancement at around 1.4 GeV 
as in the neutral B case. (We do not have enough statistics here to fit the Mb distribution 
in bins of cutt" mass.) The cutt" fraction of the (47r)~ final state is 25%, and all the uJ7i~ is 
consistent with coming from the A~ . 



VI. ANALYSIS OF D*+ujtt- DECAY ANGULAR DISTRIBUTIONS 

The is produced along with a spin-1 D* from a spin-0 B. If the A~ is spin-0 the 
D* would be fully polarized in the (J, Jz) = (1,0) state. If the A~ were to be spin-1 any 
combinations of z-components would be allowed. It is natural then to examine the helicity 
angle of the D*~^ by viewing the cosine of the helicity angle of the tt"*" with respect to the B 
in the D*~^ rest frame. 

Another decay angle that can be examined is that of the ujtc system. If the A~ is spin-0, 
the LJ is polarized in the (1,0) state and may be if the A~ is spin-1. Here the helicity angle 
is defined as the angle between the normal to the u decay plane and the direction of the 
A~ in the a; rest frame. For a spin-0 A~ the distribution will be cosine-squared. Again 
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FIG. 10. The invariant mass spectra of wvr" for the final state D*~^tt~^-k~'it~-k^ for all three 
decay modes, (left) Signal events satisfy cuts on B mass, AE and w mass (see text). The solid 
histogram is the background estimate from the AIb lower sideband and the dashed histogram is from 
the AE sidebands; both are normalized to the fitted number of background events, (right) The mass 
spectrum determined from fitting the Mb distribution and fit to a non-relativistic Breit-Wigner 
function that gives a mass of 1432ib37 MeV and a width of 376ib47 MeV. 



full polarization is possible if the A~ is other than spin-0, but any distribution other than 
cosine- squared would demonstrate that the spin is not equal to zero. 

For this analysis we use all three final states for the D*^ final state. To find the 
distributions we fit the number of events in the Mb candidate plot selected on different 
angle bins. The utt mass is required to be between 1.1 and 1.9 GeV. This restriction leaves 
111±13 events. 

In Fig. |l^ we show the helicity angle distribution, cos6d* for the D* decay. The data 
have been corrected for acceptance. We also show the expectation for spin-0 from the Monte 
Carlo. The data have been fit for the fraction of longitudinal polarization. We find 

Y = 0.63 ±0.09 . (13) 

The systematic error is much smaller than the statistical error. 

The cos 9d* distribution is not consistent with full polarization, yielding a of 17.7 for 
5 degrees of freedom. The helicity angle distribution for the A~ —>■ cuvr", cos^^;, is shown on 
Fig. |l^. Furthermore the cos 9^ distribution is quite inconsistent with a cos"^ O^j, yielding a 



X of 109 for 5 degrees of freedom. Therefore, we rule out a spin-0 assignment for the A~. 

To determine the we need a more well defined final state. This is provided by analysis 
of -B ^ DujTr~ decays. 



VII. OBSERVATION OF B ^ Dwvr DECAYS 
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FIG. 11. The B candidate mass spectra for the final state D*'^7r+7r^7r^7r'^, with K^tt^ . 

Signal events satisfy cuts on B mass, AE and to mass (see text), (a) for AE sidebands and (b) for 
AE consistent with zero. The curve in (a) is a fit to the background distribution described in the 
text, while in (b) the shape from (a) is used with the normalization allowed to fioat and a signal 
Gaussian of width 2.7 MeV is added. 

A. B candidate selection 



or D 



K 



Here we study the reactions B Dujti^ , with either a K^n^ 
decay. Other or decays have substantially larger backgrounds. 

Although we are restricting our search to cu's, we define two n~^7r'TT^ samples. One 
within 20 MeV of the known ut mass (782 MeV) and the other in either low mass or high 
mass sideband defined as three tt mass either between 732 and 752 MeV or between 812 and 
832 MeV. We also require a cut on the u Dalitz plot of r < 0.7. 

To reduce backgrounds we define 



AM, 



D 



AM,.. 



a{AM^ 



D 



a{AM^ 



AM^o 



(14) 



where AMd is the invariant candidate mass minus the known mass, AM^ is the 
invariant candidate u) mass minus the known u mass, and AM^o is the measured 77 invariant 
mass minus the known 7r° mass. The a's are the measurement errors. We select candidate 
events requiring that xl is < 12 for the Kir mode and <6 for the Ktctt mode. 



B. B- D^uTT- Signal 

We start with the D° K^tt^ decay mode, for events in the u peak. We show the 
candidate B mass distribution, Mb, for AE in the side-bands from -7.0 to -3.0a and 7.0 
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FIG. 12. The invariant mass spectra of vr+vr^vr'^ for the final state D*^tt^tt^tt^tt^ for the 
K^TT^ decay mode. The solid histogram is the background estimate from the Mb lower 
sideband and the dashed histogram is from the IS.E sidebands; both are normalized to the fitted 
number of background events. There are two combinations per event. 



to 3.0cr on Fig. This gives a good representation of the background in the signal 

region. The /S.E resolution is 18 MeV (cr). We fit this distribution with the shape given in 
Equation |. 

We next view the Mb distribution for events having Ai? within 2a around zero in 
Fig. |T^(b). This distribution is fit with a Gaussian signal function of width 2.7 MeV and the 
background function found above whose normalization is allowed to vary. We find 88 ±14 
events in the signal peak. 

We repeat this procedure for events in the u sidebands. We use for our definition 
pseudo-cij masses in the center of the sideband intervals. We show the Mb distribution for 
events in the AE sideband, defined above, and those having AE within 2a around zero in 



Fig. 18. We find no significant signal. 



D^wK Signal 



The same procedure followed for the final state is used for the final state. We 
show the candidate B mass distribution, M^, for events in the AE side-band on Fig. |19|(a). 
The AE resolution is 18 MeV [a). This gives a good representation of the background in 
the signal region. We fit this distribution with a shape given in equation ^. 

We next view the Mb distribution for events having AE within 2a around zero in 
Fig. |l^(b). This distribution is fit with a Gaussian signal function of width 2.7 MeV and the 
background function found above whose normalization is allowed to vary. We find 91±18 
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FIG. 13. The Mb spectra for D*^uj'k~ for the K~tt~^ decay mode, (a) AE sidebands 
and (b) AE around zero. 



events in the signal peak. 

We repeat this procedure for events in the u sidebands. We show the Mb distribution 



for both AE sidebands and AE within 2a around zero in Fig. |20 



There is no evidence of any signal in the lj sideband plot, leading to the conclusion that 
the signal is associated purely with u. 



D. Branching Fractions 

We determine the branching ratios, shown in Table by performing a Monte Carlo 
simulation of the efficiencies in the two modes. We use the current particle data group 
values for the relevant u, and branching ratios of (88.8±0.7)% [uj — >■ 7r"'"7r~7r°), 
(9.0±0.6)% iD+ K-n+n+) and (3.85±0.09)% (D^ K-ti+) [g. The efficiencies listed 
in the table do not include these branching ratios 0]. 



TABLE III. Branching Fractions for the Dujtt final state 



D Decay Mode 


Fitted 7^ of events 


Efficiency 


Branching Fraction (%) 




88±14 


0.064 


0.41±0.07±0.06 




91±18 


0.046 


0.28±0.05±0.04 



The systematic error arises mainly from our lack of knowledge about the tracking and tt*^ 
efficiencies. We assign errors of ±2.2% on the efficiency of each charged track, and ±5.4% 
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FIG. 14. The invariant mass spectra of ujtt^ for the final state D""tt ' it tt it 

K^TT^ decay mode, (left) Signal events satisfy cuts on B mass, AE and u mass (see 
text). The solid histogram is the background estimate from the Mb lower sideband and the dashed 
histogram is from the AE sidebands; both are normalized to the fitted number of background 
events, (right) The data fit to a non-relativistic Breit-Wigner signal and a smooth background 
function. The mass and width are 1367ib75 MeV and 439ibl35 MeV, respectively. 



for the 7r°. The error due to the background shape is evaluated in three ways. First of all, 
we change the background shape by varying the fitted parameters by la. This results in 
a change of ±5.0%. Secondly, we allow the shape, p2, to vary (the normalization, pi, was 
already allowed to vary). This results in 5.5% increase in the number of events. Finally, we 
choose a different background function given in equation ^ and repeat the fitting procedure. 
This results in a 1.0% decrease in the number of events. Taking a conservative estimate of 
the systematic error due to the background shape we arrive at ±5.5%. 



E. The ujTT System 

For all subsequent discussions we add the and final states together. We select 
sample of cj's in the 7r"'"7r~7r° mass window of 782±20 MeV using only combinations having 
r < 0.7 in the Dalitz plot. 



In Fig. ^we show the un~ mass spectrum in the left-side plot. The solid histogram shows 
events from the lower Mb sideband region (5.203 - 5.257 GeV) suitably normalized. The 
dotted histogram shows the background estimate from the AE sidebands, again normalized. 
In the signal distribution there is a wide structure around 1.4 GeV, that is inconsistent with 
background. We re-determine the 0J7r~ mass distribution by fitting the Mb distribution in 
bins of uJ7r~ mass, and this is shown on the right-side. 

This structure appears identical to the one we observed in B ^ D*tun^ decays. 
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FIG. 15. The cosine of the angle between the and the D* flight direction in the D* rest 
frame for the D*A~ final state (solid points) along with a fit (solid curve) allowing the amount of 
longitudinal and traverse polarization to vary. The dotted curve is the expectation for a spin-0 A~ . 



F. Angular Distributions in Dlutt 

We can determine the spin and parity of the A~ particle by studying the angular dis- 
tributions characterizing its decay products. The decay chain that we are considering is 
B ^ A D; A UTT and u tt+tt^tt^. The helicity formalism [|ll|] is generally used in the 



analysis of these sequential decays. This formalism is well suited to relativistic problems 
involving particles with spin J and momentum p because the helicity operator h = J ■ p is 
invariant under both rotations and boosts along p. 

There are two relevant reference frames. The first one, that we will define xai/aZa is the 
rest frame of the A particle, with the axis pointing in the A direction of motion in the B 
rest frame. The x'a direction is arbitrary. The second one, x^^y^Zuj, is related to xaVa^a by 



the rotation through 3 Euler angles (f)A,dA, —4>a, as shown in Fig. The angle 0^ defines 
the orientation of the plane containing the u direction in the A rest frame and the 2C1 axis 
with respect to the x^a — za plane. The angle 6a is the polar angle of the u momentum 
vector in the A rest frame. Note that the A decay plane has an azimuthal angle (pA both in 
the xaUaZa and in the x^^y^z^^ references. The angles 9^^ and define the orientation of the 
(jU decay plane in the u rest frame. As the angle (pA is arbitrary, the only angle that has a 
physical meaning is x = — 0a;, the opening angle between the A decay plane and the u 
decay plane. 

Both the B meson and the D meson are pseudoscalar, therefore their helicity is 0. Thus 
A will be longitudinally polarized independently of its spin. In order to calculate the decay 
amplitude for this A uj-k^ process, we need to sum over the uj helicity states: 



here dH^J 



A = J:xMt('pA,eA,-^A)Dlloi<P^,9^,-<P^)B,^o, (15) 
—4>a) is the rotation matrix that relates the xaVaza and the x^y^z^ frames 
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FIG. 16. The cosine of the angle between the normal to the tt'^tt^tt decay plane and the lu 
boost direction for the D*A~ final state (solid points) along with a fit (solid curve) allowing the 
amount of longitudinal and traverse polarization to vary. The dotted curve is the expectation for 
a spin-0 A~ . 

and D'^l^Q^cpi^jO^, —(pi^) is the rotation matrix relating the Xi^^y^Zi^ and the direction of the 
normal to the u decay plane n{6^, 0^). 

In general, there are three helicity amplitudes that contribute to this decay: Biq and 
-B_io, corresponding to a transverse u polarization, and Bqq, corresponding to a longitudinal 
uj polarization. This expression can be simplified by observing that A ujti is a strong 
decay and thus conserves parity. Thus, the helicity amplitudes are related as: 



where S{A) is the spin of particle A and rjA, rj^ and t/t, represent the intrinsic parity of the 
decaying particle and its decay products, respectively. 

Eq. |l^ relates the two transverse helicity amplitudes, while Eq. |1^ forbids the presence 
of a longitudinal component under certain conditions. For example, if A is a 1~ object, the 
uj has transverse polarization and -B_io = —Biq. When the sign in Eqs. |16|-0 is positive, 
two parameters determined by the hadronic matrix element affect the angular distribution 
and thus we cannot fully determine it only on the basis of our assumptions on the A spin 
parity. We have carried out the calculation of the predicted angular distributions including 
spin assignment for A up to 2. The predicted angular distributions are summarized in Table 



The statistical accuracy of our data sample is not sufficient to do a simultaneous fit of 
the joint angular distributions shown above. Thus only the projections along the 6a, Ow and 




(16) 




(17) 



rvi. 
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TABLE IV. Differential angular distributions (modulo a proportionality constant) predicted 
for different spin assignments. (Note 0"*" is forbidden by parity conservation.) 



da /d cos 6 Ad cos 6 i^dx 
0^ l-Boopcos^^^ 
1 ~ I -Sio P sin^ 6a s\t? sin^ x 

1+ l-^ioP sin^ 6a sin^ duj"^ cos + I-SqoP cos^ 6a cos^ d^j 

- 1 /2Re{BioB^Q) sin 26 a sin 29^ cos x 

2- 3 1 5io P sin2 261a sin^ 6^ cos^ x + | Bqo P (3 cos^ - 1 ) ^ cos^ 6, 

- VSReiBioB^o) sin 26'a(3 cos^ 9a - 1) sin 29ijj cos x 
2+ 3/4 1 Sio P sin^ 29a sin^ 0^ sin^ x 



TABLE V. Projection of the angular distributions along the cos^a, 

cos 9^^ and x axes. 



J da/dcos9A da/dcos9u) da/dx 

0= f |SooP 47r|Boopcos2 0^ 4/3|Soop 

1- f iSiopsin^^A f l^iopsin^^^ fl^iopsin^x 
1+ f (iSiopsin^^A ^{\B,o\^sm^6^ |(4|i3io|2 cos^ x 

+ |SoopCOs2 0A) +|SoopCOs2^^) +|5oo|') 

2- f (3|5iopsin2 2eA if^d^iopsin^e^ 4|Sio|2cos2x 

+|SooP(3cos2 0a-1)^) +\Boo\^ cos^ 6^) +\Boo\^ 

2+ 7r|Biopsin2 20A flBiopsin^^^ IflBiopsin^x 
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FIG. 17. The B candidate mass spectra for the final state D^ujtt~ , with K~Tr~^. (a) for 

AE sidebands, and (b) for AE consistent with zero. The vertical scale in (a) was multiplied by 0.5 
to facilitate comparison. The curve in (a) is a fit to the background distribution described in the 
text, while in (b) the shape from (a) is used with the normalization allowed to float and a signal 
Gaussian of width 2.7 MeV is added. 



X are fitted, integrating out the remaining degrees of freedom. Table |V], gives the analytical 
form for these projections. 

We determine the projections of these angular distributions by fitting the Mb distribution 
as a function of the various angular quantities cos^^, cos6'^, x- We restrict the U7i~ mass 
range to be between 1.1 and 1.7 GeV, containing 104 signal events. In order to fit the 
angular distribution with theoretical expectations, we must correct the data for acceptances. 
We determine the acceptance correction by comparing the Monte Carlo generated angular 
distributions with the reconstructed distributions. The angular dependent efficiencies are 



shown in Fig. ^ 



The corrected angular distributions are shown in Fig. The data are fit to the expec- 
tations for the various assignments. For the 0~, 1~ and assignments, the curves have 
a fixed shape. For the 1"*" and 2~ assignments we let the ratio between the longitudinal and 
transverse amplitudes vary to best fit the data. We notice that the uj polarization is very 
clearly transverse (sin^^^) and that infers a 1^ or 2^ assignment. 

We list in Table |V| the x^ /dof for the different assignments. The 1~ assignment is 
preferred, having a x^/dof of 1.7. The other assignments are clearly ruled out. If the 1~ 
assignment is correct, the probability that our fit yields a x^ jdof equal to 1.7 or greater is 



3.8% p 
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FIG. 18. The B candidate mass spectra for the final state D^ujtt~, with K^ir^ and lu 

sidebands for AE sidebands (histogram) and AE consistent with zero (points). The AE sideband 
numbers have been divided by 2. 



TABLE VI. Results of fits to angular distributions 





0- 


1+ 


1" 


2+ 


2- 


xVdof 


7.0 


4.5 


1.7 


3.2 


5.3 


dof 


15 


14 


15 


15 


14 


probability 


1.9 X 10-15 


3.3 X 10-s 


3.8% 


2.7 X 10-5 


3.3 X IQ-io 



VIII. DISCUSSION OF NATURE OF THE A' 

We have found a 1" object decaying into ujtt'. A non-relativistic Breit-Wigner fit assum- 
ing a single resonance and no background gives a mass around 1420 MeV with an intrinsic 
width about 400 MeV. Signals for un^ resonances have been detected before below 1500 
MeV. There is a well established axial- vector state, the bi(1235), with mass 1230 MeV and 
width 142 MeV. Data on vector states, excited p's, are inconsistent. Clegg and Donnachie 
13| have reviewed —>■ {A7c)~u, e~^e~ — > h'^tt' and e^e~ — >■ Tr+Tr+vr-Tr- data, including 



the un final state. Their best explanation is that of two 1 states at 1463±25 MeV and 
1730±30 MeV with widths 311±62 and 400±100 MeV, respectively. Only the lighter one 
decays into ujtt. The situation is quite complex, however. They conclude that these states 
must be mixed with non-gg states in order to explain their decay widths. There is also an 
observation of a wide, 300 MeV, cutt" state in photoproduction at 1250 MeV [T^, that is 



dominantly the 6i(1235) with possibly some 1 in addition. Our state is consistent with 
the lower mass p'. We do not seem to be seeing significant production of the higher mass 



24 



1630301-010 



100 
75 
50 



« 25 



S> 100 



75 
50 
25 



(a) 




5.20 



5.22 



5.24 



5.26 



5.28 



5.30 



FIG. 19. The B candidate mass spectra for the final state D~^ujtt~, with K^ir^ir^ (a) 

for AE sidebands and (b) for AE consistent with zero. The vertical scale in (a) was multiplied 
by 0.5 to facilitate comparison. The curve in (a) is a fit to the background distribution described 
in the text, while in (b) the shape from (a) is used with the normalization allowed to float and a 
signal Gaussian of width 2.7 MeV is added. 



state into ujn~, as expected. 

Several models predict the mass and decay widths of excited p and u mesons . For 
example, according to Godfrey and Isgur [|I6| the first radial excitation of the p is at 1450 
MeV. There is a large variation among the models, however, on prediction of the relative 



decays widths ranging from no tttt to mr being equal to unr |]T7|. 

Since we have observed a wide 1^ state in the mass region where the p' is expected, the 
most natural explanation is that we are observing the p' for the first time in B decays. 

We note that r~ lepton decays into un~' have been observed, and the 1~ spin-parity 
definitely established. However, the relatively low mass of the r~ distorts the the mass 



spectrum significantly, and makes it difficult to extract the p' mass and width |18| |19 . 



IX. MASS AND WIDTH VALUES FOR THE p' 

Here we find the best values for the p' mass and width. This procedure is discussed in 
more detail in the Appendix. The shape of the uj7r~ mass spectrum is affected by the phase 
space allowed by the B decay, the B decay amplitude, the decay amplitude for the p', and 
finally the shape of the Breit-Wigner decay distribution. For B D^*^ujn^ decay: 

dT{B dWcjtt) = -J— \A{B D^*^p') BW{p') A{p' cj7r)p (18) 
X dV{B Z}(*V') dV{p' ^ uj-n) , 
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FIG. 20. The B candidate mass spectra for the final state D~^ujTr~, with K~7r^ir~^ 
and uj sidebands for AE sidebands (histogram) and AE consistent with zero (points). The AE 
sideband numbers have been divided by 2. 



where D^*'' indicates either a D* or a D meson. 

The phase space for two-body decays is well known. The decay amplitude for the B 
decay can be obtained from factorization where the p' is assumed to be identical to the 
lepton current [Q. Finally the Lorentz structure of the p' decay can be accounted for and 
we are left only to consider a Breit-Wigner amplitude of the form 

Breit - Wig.e.,M„.) ^ _ g^^'^lp,^,) ■ 
where the mass dependent width is given by 

r(M )-r(M) (pAM^^X f^p'Y 1 + iRpUM,')Y 
T{M^.) - r(M,) [jj^j (^^j 1 + {RpUm^^W • 

We allow the parameters R and n to float in the fit. They represent parameterizations of the 
hadronic matrix element and Blatt-Weiskopf damping factors. R is the p' radius in units of 
fm/hc. Fig. ^ shows the fit to M^t^ distribution, where we have summed the D^, and 
D*^ data. (We have corrected the data in each channel for the M^jr efficiency dependence, 
which is small for the and modes.) 

The fit gives the mass and width to be (1349±25l^°) MeV and (547d=86lt^) MeV, re- 
spectively. The errors are statistical and systematic; they are derived in the Appendix. The 
values for R and n are 2.05lj';i| fm/^c and 0.57j;o;83, respectively. 
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FIG. 21. The invariant mass spectra of wvr" for the final state DLijn~ for both D decay 
modes, (left) The solid histogram is the background estimate from the Mb lower sideband and 
the dashed histogram is from the AE sidebands; both are normalized to the fitted number of 
background events, (right) The mass spectrum determined from fitting the Mb distribution and fit 
to a Breit-Wigner function. We find a peak value of 1415±43 MeV and a width of 419±110 MeV. 



X. SEARCH FOR OTHER RESONANT SUBSTRUCTURE IN D*(47r)- 

We have accounted for ~20% of the (47r)~ final state. We would like to disentangle other 
resonant substructure. Since the background is large in modes other than —>■ K'tt^ we 
will only use this mode. One process that comes to mind is that where the virtual W~ 
materializes as an a^", that subsequently decays into tt+tt^tt" and we produce a that 
decays into D*~^7r^. This process should be the similar to that previously seen in the reaction 



B- 



D**^n , where the D**^ decayed into a D 



TX 



We search for the presence of an 



aj~ by examining the ti^ti^ti^ mass spectrum in Fig. 

There is an excess of signal events above background in the aj~ mass region, that cannot 
be definitely associated with the ai. Proceeding by selecting events with tt^tt~tt~ masses 
between 0.6 and 1.6 GeV, we show the D*"'"7r° invariant mass spectrum in Fig. p7 . 



Although there is a suggestion of a low mass enhancement, it is not consistent with D** 
production that would peak in region of 2.42 - 2.46 GeV. Perhaps we are seeing an indication 
of fragmentation at the b c decay vertex here. 

We also display for completeness the "aj~7r°" mass distribution in Fig. There may or 
may not be a wide structure in the (An) 
substructure in this decay channel. 



mass. At this point we abandon our search for 
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FIG. 22. Relationship between the A rest frame xaVaza and the uj rest frame x^y^z^. xa and 
x^ he in the same plane. 
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FIG. 23. Reconstruction efficiency dependence on (a) cos^^i, (b) cos^^^, and (c) x- 



XI. DISCUSSION AND CONCLUSIONS 



We have made the first statistically significant observations of six hadronic B decays 



shown in Table VII 



There is a low-mass resonant substructure in the ujti mass. The fit to a sophisticated 
Breit-Wimer function eives the mass and width to be (1349±25l^°) MeV and (547±86l^^) 
MeV, respectively. 

The structure at 1349 MeV has a spin-parity consistent with 1^. It is likely to be the 
elusive p' resonance [l^. These are by far the most accurate and least model dependent 
measurements of the p' parameters. The p' dominates the final state. (Thus the branching 
ratios for the D^*'^ujtx~ apply also for D^*^p'~.) 

Heavy quark symmetry predicts equal partial widths for D*p' and Dp'. We measure the 
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FIG. 24. The angular distribution of 9a (top-left), 9^^ (top-right) and x (bottom). The 
curves show the best fits to the data for different assignments. (The 0~ and l"*" are almost 
indistinguishable in cos 9a, while the 1~ and 2+ are indistinguishable in cos9u and x-) 
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FIG. 25. Fit to the M^^ distribution. The data in the D , D'^ and D*^ channels have been 
summed and corrected using mass dependent efficiencies. 



TABLE VII. Measured Branching Ratios 



Mode 


B (%) 


# of events 




1.72±0.14±0.24 


1230±70 


S° ^ D*+LJ7r- 


0.29±0.03±0.04 


136±15 


S° ^ D+LOTT- 


0.28±0.05±0.04 


91±18 


B- D*^-K+TT-TT-TT^ 


1.80±0.24±0.27 


195±26 


B- D*°LJ-K- 


0.45±0.10±0.07 


26±6 


B' D"UJTT- 


0.41±0.07±0.06 


88±14 



r( 




^ D* 


V'-) 


r 









relative rates to be 

r f# ^ D*+o'-'] 

1.04 ±0.21 ±0.06 (21) 



-\ = 1.10 ±0.31 ±0.06 (22) 

1.06 ±0.17 ±0.04 . (23) 



T[B ^ Dp'- 

Thus the prediction of heavy quark symmetry is satisfied within our errors. 

Factorization predicts that the fraction of longitudinal polarization of the D*^ is the 
same as in the related semileptonic decay B — > D*i~u at four-momentum transfer equal 
to the mass-squared of the p' 
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FIG. 26. The invariant mass spectra of tt+tt tt for the final state Z?*+7r+7r vr tt^ for 
D° K--K+. The solid histo gram is the background estimate from the Mb lower sideband 
and the dashed histogram is from the AE sidebands; both are normalized to the fitted number of 
background events. 



t(b ^ D*+p'-) V (b 



D* 



V 



9 2 

p 



(24) 



Our measurement of the D*^ polarization (see Fig. 15) is (63±9)%. The model predic- 
tions in semileptonic decays for a of 2 GeV^, are between 66.9 and 72.6% I^Tj. Thus this 
prediction of factorization is satisfied. 

We can use factorization to estimate the product of the p' weak decay constant fpi and 
the branching ratio for p'~ — > cuvr". The relevant expression is 



v(b^d*^p'-, p'- 



LUTT 



D* 



Q^^clfl,B ip' 



IK,, 



'ud\ 



(25) 



where Ci is a QCD correction factor. We use ci=l.l±0.1 [p^ . 



We use the semileptonic decay rates given in Barish et al. [23]. The product 

fpS [p - con-) = 0.011 ± 0.003 GeV^ , (26) 



where the error is the quadrature of the experimental errors on the experimental branching 
ratios and ci. 

The model of Godfrey and Isgur predicts decay constants widths and partial widths of 
mesons comprised of light quarks by using a relativistic treatment in the context of QCD 
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FIG. 27. The invariant mass spectra of L>*+7r for vr+vr vr masses between 0.6 - 1.6 GeV 



for the final state D* 



TT TT TT vr 



with D 







K 7r+. The soHd histogram is the background 



estimate from the Mb lower sideband and the dashed histogram is from the A£^ sidebands; both 
are normahzed to the fitted number of background events. 



Tq] . They predict both fpt and B {p' 



The branching ratio prediction is beheved to be more accurate |24 



Lun ); the values are 80 MeV and 39%, respectively. 

We use this to extract 



fp> = 167 ± 23 MeV 



(27) 



The model predicts a lower value for fpi than observed here, if factorization is correct. 

We note that all the B — > D^*^p' branching ratios that we have measured are approxi- 
mately equal to the B D^*^p branching rates @] if a model value of B (p'^ un^) = 39% 
is used. 

Finally, although the B — * D*+(47r)~ and B~ D*°(47r)~ branching ratios are nearly 
equal, the cuvr" branching ratios are about 1.5 times larger for the charged B than the neutral 
B, maintaining the trend seen for the tt~ and p~ final states. Since the B~ lifetime is if 
anything longer than the B^, this trend must reverse for some final states. It has not for 
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FIG. 28. The invariant mass spectra of Tr+vr^vr^vr for Tr+vr^Tr^ masses between 0.6 - 1.6 
GeV for the final state D*^TT^TT~7r~7r^ with K~'k^ . The solid histogram is the background 

estimate from the Mb lower sideband and the dashed histogram is from the A£^ sidebands; both 
are normalized to the fitted number of background events. 
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XIII. APPENDIX: DETERMINATION OF THE p' MASS AND WIDTH 



A. Introduction 



Since the p' decays to tuvr" via a P-wave, we need to take into account the fact that 
the width, r(M^^), may not be constant, but can vary with M^^r- Furthermore we need to 
consider the kinematic limits from B D^*^p' decay and p' — > cutt decay. 



B. The Differential Decay Distribution 



We can write a general expression for the differential distribution for M^^t^ in dV{B D^*^uj7i) 
via the p' as: 

dT{B dWc^tt) = \A{B D^*^p') BW{p') A{p' un)]^ (28) 

X dV{B D^*^p') dV{p' uo-k) , 

2tx 

where dV indicates a phase space term, BW indicates some form of a Breit-Wigner shape 
function, and A indicates an amplitude. 

The phase space oi B ^ D^*^ p' provides a cut off at higher Ai^^, whereas the phase 
space of p' — > ojTi provides a cut off at lower M^^^,. We assume that the two decay stages are 
independent and can be factorized. Thus the function can be calculated from the widths of 
B and p' decays. 



1. The B Decay Width 



The width for p' production is given by: 
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2M 



B 



\A{B DWp')P d'PiB D*p') 



dV{B D^*^p') = — 



1 f2pD* 



(29) 
(30) 



Svr V Mb 

A{B^D'^*^p') ~ ^G^Vcfc X Lorentz structure X /p/(A4,^) 

where fp'{M^T,) is the p' weak production form factor. 

Instead of having to calculate fp/{M^T^), we can use our knowledge of semileptonic b 
decays, B — > D^*H^T7 coupled with factorization to approximate A(B — > D*p'). 



Factorization tells us: 

T{B ^ Z^Wp') = Stt^ iKdl' fJM^ 



\ai\ 



(31) 



« PdW Ml X {|i/+(Ml)P + \H.{Ml^)\' + |i/o(Ml)|n 



The helicity amplitudes H{M^^) can be related to the axial- vector form factors Ai{q ) and 
A2{q'^), and vector form factor V{q'^). Details can be found in reference p5[ . 



{Ms + M^O Ai(g^) T V{q' 



1 



2Md*V^ 



Mb + Mb* 
(M| - Ml,, - q^){MB + M^O Ai(g2) 



(32) 



^2(g') 



Mb + Mz)* 

In the heavy-quark symmetry limit, the form factors Ai, A2, and V are related to the 
Isgur-Wise function. With correction due to finite heavy quark mass and a^, they can be 
written as 



1 - 
R2 

Ri 



{Mb + Md'Y 
Mb + Md- 



Mb + Mp* 

2yJMBMD, 



(33) 



2VMbMd* 
Mb + Md- 



2^MbMd 

were w is the invariant four-velocity transfer. 



The values calculated by Neubert for Ri and R2 have the explicit dependence on w [26 



of 



1.35 - 0.22(w - 1) + 0.09(w - 1) 



(34) 



(35) 



Ri{w) 

R^(w) = 0.79 + 0.15(«; - 1) - 0.04(«; - 

while Close and Wambach determine 

Ri{w) = 1.15 - 0.07(?i; - 1) + C(w - 1 
R^(w) = 0.91 + 0.04(w-l) + C(w- 1)2 . 

CLEO measured Ri{0) and i?2(0) to be 1.18 ± 0.30 ± 0.12 and 0.71 ± 0.22 ± 0.07, respec- 
tively The form factor Ha^ (w) can be assumed to be linear as Ha^ (w) = 1 — p\_^ {w — 1); 
CLEO measured p% to be 0.91 ± 0.15 ± 0.06. 

r{B D* p') as function of M^tt is shown in Fig. where the form factor fpi is not 
included. 
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FIG. 29. Decay width of i? ^ D* p' as function of M^-,^. The parameters i?i and i?2 are from 
Neubert (sohd line), Close and Wambach (dashed Hne) and the CLEO measurement (dotted hne), 
with as 0.91. 



2. The Width of p' decay 



The width of p' decay can be expanded as: 

1 



2M. 



\A{p' ^ujTi)Y dV{p' ^cutt) 



dV\p — > ujtt) — — ' 



Lorentz structure x h{M^^) , 



(36) 

(37) 
(38) 



where is evaluated in the un rest frame, and h{M^^) is the p' strong decay form factor, 
which is usually assumed to be a constant. 

The Lorentz structure of p' ujh must be linear in the polarization vectors of the p' 
and the uj. The simplest mathematical expression is 6p> ■ St^, which is S-wave. For P-wave, 
since Epi ■ Pp' — and e^, ■ Poj — hy transversality, the only Lorentz scalers we can form are 
{sp/ ■pu)){£u) • Pp') and e^,/a/j ep,e';^p^Pp,. The first term violates parity conservation. Later, we 
describe the detailed calculation of the second term, which gives 



Thus, the width of p' is: 



r(p'^uj7r) (X h'^{M^ 



(39) 



(40) 



Note, that the contribution of longitudinal polarized uo to the width is zero, which is 
consistent with our spin-parity study; our measurements gave TL^g/Ttot to be (10 ± 9)% in 
the £)*a;7r mode and (—0.4 ± 22%) in the Duj-k mode. 
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3. The Breit-Wigner 



The Breit-Wigner function has a long history. Fundamentally we are approximating the 
decay as having a non-changing amplitude as a function of mass in the simplest case. In 
general the denominator of the Breit-Wigner has a fixed form in amplitude given by 

Denominator = M^^ - - iM^^Ttot{M^^) , (41) 

The numerator is where we have to have more discussion. Since uott is a large decay 
mode of p', we assume that the mass dependence of Vtot can be approximated by the mass 
dependence of F^^ 

TUM^^)=TUM,)^-f^^^ (42) 



From the derivation of F(p' un) above, we see that 

r^.(MpO 




(43) 



4- The Differential Distribution 
Eq. ^ now can be rewritten as: 

.r(B ^W, . TiB ^V) „^^_g;;^-gji(,^., 1g- (44) 

With the relation shown in Eq. we have the differential function which can be used 
for fitting the mass distribution: 

where F is the M^^-dependent total width. From Eq. ^ and Eq. ^ we get that: 

^ = ^„xfM^elyxf^^^^V (46) 



5. Decay Form Factor of p' 



uj-n 



Now, what form should we use for h{M^^)7 In Eq. the dimensions of h{M^^) are 

GeV""*^. But the dependence on M^^vr could be anything. So we can try h{M^^) oc M"^, 
where n is allowed to float in the fit. 

h{M^^) could also include Blatt-Weisskopf factors p9|. A Blatt-Weisskopf factor for the 
P-wave decay p' ljti would be of the form: 

where R is the radius of the p' meson. A typical value is = 1 im/hc. 
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C. Fit to M^TT Spectrum 

Since we are limited by the statistics, we chose to add the D and D* final states together. 
This is permissible since the Lorentz structure of p' decay in the different modes is the same. 
The only difference is in the width oi B ^ D^*^p' part, that depends on M^^, however this 
difference is slight. 

We use following formula to represent the width: 

r(M„,, ^ r(M,) f^y" ^^^l^4M. (48) 



1 + (Rp^M^,))- 

We allow the parameters R and n to float. We use Neubert's calculated values for Ri and 



R2; effects of changing this to other estimates will contribute to the systematic error. Fig. |25 
shows the flt to M^^^ distribution. 

The flt gives the mass and width to be (1349±25) MeV and (547±86) MeV, respectively. 
The values for R and n are 2.05l[j4 ^md 0.5711qII, respectively. 

Systematic errors can arise from several sources. We test our flt globally by restricting 
the mass range to below 1.8 GeV. This results in a shift in the mass by -5.7 MeV and the 
width by -23.2 MeV. Both of these changes are much smaller than the statistical error. 

Two speciflc sources of systematic error are using the D* mass to represent both D* and 
D meson flnal states, and using the form-factors -Ri and R2 from Neubert's calculation. To 
estimate the error on the former, we use the D in the flt; for the latter we use form-factors as 
calculated by Close and Wambach or, the CLEO measurement of these parameters. Listed 
in the Table |V111| are the contributions from each of these sources. 

Finally, to estimate the error due to our model of the cuvr" shape, we allow the values 
of R and n to change so that the of the flt increases by one unit, corresponding to a 



one standard deviation variation. We list in Table |VIII| the maximum positive and negative 



changes in the mass and width values allowed by these variations. Thus our flnal value for 
the mass of p' resonance is (1349 ± 25l5°) MeV. The width is (547 ± 86l4^) MeV. 



TABLE Vlll. Mass and width fit of resonance p' 


Parameters 


Mass (MeV) 


Width (MeV) 


Neubert (nominal) 
M = Md 
Close- Wambach 
CLEO D*lu 
Vary R and n 


1349 ± 25 
-4.0 
-1.5 

+8.2 

+5.8 
-3.2 


547 ± 86 
-31 
-9 
+24 

+40 
-31 


Systematic error 


+10.0 
-5.3 


+46 

-45 



It is possible that values of the p' mass and width could be affected if additional resonance 
substructure of unknown origin were to be included. There is no evidence, however, that 
such structures are needed. By including all the known physics effects of phase space and 
Lorentz structure we are able to describe the data quite well. 
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D. Calculation of the Lorentz Structure 



In this Section, we will calculate the Lorentz structure le^^^a/s ^p'^J^P^Q'^P) where we use 
and to represent the four momentum of p' and u. In general the helicity structure is 

more apparent when calculated in a rest frame where the p' is in motion. Later we will give 

the expression in p' rest frame. 

The transversality ■ = |^. For spin 1 massive particle with momentum p = pcz, 

the helicity states can be represented as: 

£(±) = T(0,l,±z,0) (49) 
= {p,0,0,E)/M 

For general p, we can construct e(^o) as follows. The e(^±) are not shown as they will not 
be directly used. 

- (irr If) (-) 

Now we need to prepare calculations of several quantities. First is gp^-^gp^-,. The summa- 
tion over all helicity states is given in PD[. For helicity 0, we use Eq: and expand it. We 



have 



E = -9^' + ^ (51) 



A=o,± 



+ = -^7^^^ + 5.0^.0 - (53) 

In {pf^p")*, only /i, i/ = 1,2,3 have non-zero values. In the later calculations, we always 
assume that transverse helicity states (A = itl) have the same probabilities. 

For terms like eaji-ysp^p^ ■ ■ ■, exchanging a and /? will introduce a negative sign. Thus the 
value is zero. 

ea/s-ys p'^p^ X anything = (54) 

The Lorentz structure will be calculated with summation of all helicity states and different 
polarizations. 

Atot = Aqo + Aoi + Aio + Au = e«/375 ^kX^^u (^"^"pV)^' {e'^e^q^p")^ (55) 

where Aoi for longitudinal polarized p' and transverse polarized u. Others are similar. 
There is one term which often shows up; expanding we find: 

eoijkeoimn PjPm qkq-n = plql + plql - 2p2P3g2g3 {i = 1) (56) 

+ pIqI + phi - "^pmqiqs {i = 2) 
+ pIqI + pWi - 2piP2gig2 {i = 3) 
+ {plql + plql + phD - (plql + pIqI + pIqI) 

= ipl+Pl+Pl){(ll + (l2 + (ll) - {PiQi + P2q2 + PsQaY 
= p^(f — {p- qf' = p^q^{'i — cos^O) 
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Now the summation over all helicity states: 

Atot = e^p,s e.A^. + ^) pV (-9^' + q'q" (57) 



M2 



= -2eo/3^5 eoA/xr. p'^p^q^q" {a = or k = 0) 
+ eao^is eaOf,u p^p^q^q" (/3 = A = 0) 
+ eaiBuo (^aXfiO J^p^q^q^ (5 = zy = 0) 
+ 2eaOM5 e^AMO p^p^q^q"" (/? = i/ = or A = 5 = 0) 

= -2p2g2(l - cos'^9) + 2Ejg2 + 2^^^^ - AEp>E^pq cosO 

Lorentz structures for longitudinal polarized p' and u: 

Aoo = 6.,,,6.A.. - 5.0^.0 + /P^ - 5.0^.0 + q'q^ (58) 

= ea/375e.AM. {pyyp^P^ (q"fq>^yqSq'^J_ 



p2q2 



= 

In the first step, only one of nine terms survive; seven terms are eliminated due to Eq. 
The one which has four 5's results in two subscripts in each e. In the second step, Eq. ^ 
limits that f3, A, 6 and u are all zero, which also gives zero. 

The Lorentz structures for longitudinal polarized p' and transverse polarized u are: 

Aoi = e^f,^se.x,u - 5.o5«o + pV (^-9^" + '^^o'^mo - q'q" (59) 

= eap-ys e^A/ii. f^aof^KO 9^^ p'^p^q^q" - eaf3ys e^x^iu 9^^ — — p'^p^q^q" 

— —^0|3^lS ^OXfiu P P q q + ^aOfiS ^nOfiu — — P p q q 

= - 1) pY{1 - cos'9) 

= M}q\l - cos^e) . 

Similarly, Lorentz structures for transverse polarized p' and longitudinal polarized uj are: 

^10 = mlp\l-cos^e) , (60) 

while Lorentz structures for transverse polarized p' and uj are given by: 

All = Atot — Aqq — Aqi — Aio 

= 2Ejg2 + 2E2p2 _ AEp.E^pq cos9 - {2p'^q^ + M^q^ + mlp'^){l - cos'^O) 

= {El + EIp") (1 + cos^d) - AEp,E^pq cosO . (61) 

In the rest frame of p', the Lorentz structures are expressed as: 
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Aot = -iM^ql (62) 
Ao = 

= M^2,g2 (1 _ cos^^) 

Aio = 

An = Mlqlil + cosH) . 
(Note, that there is no contribution from longitudinal polarized a;.) 
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